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^  ABSTRACT  ■ — ' - L_ 

A  theory  is  developed  that  describes  the  processing 
of  data  collected  with  a  satellite-borne  or  airborne 
synthetic-aperture  radar  (SAR).  A  description  of  the 
target-radar  geometry,  and  the  form  of  the  received, 
demodulated  radar  signal  is  given.  It  is  shown  that 
the  solution  for  the  case  in  which  the  antenna  is 
directed  perpendicular  to  track  (sidelooking),  is 
obtained  by  a  simplification  of  the  general  case  in 
which  the  antenna  is  squinted  with  respect  to  the 
perpendicular  to  track. 

A  mathematical  description  of  the  signal  proces¬ 
sing  operations  required  to  produce  a  SAR  image  from 
the  received  radar  signal  is  presented  and  the  form 
of  the  processed  signal  is  described.  In  particular, 
a  technique  which  employs  two-dimensional  matched 
filtering  to  produce  the  radar  image  is  discussed  and 
the  ability  of  this  approach  to  accommodate  the 
coupling  of  the  range  (across-track)  and  azimuth 
(along-track)  signals  is  investigated.  In  addition, 
the  extensions  to  the  theory  required  for  non-coherent 
averaging  are  included. 

/v 

1.  INTRODUCTION 

In  this  report  a  description  of  the  mathematical  operations  required 
to  produce  geometrically  correct  images  from  airborne  or  satellite-borne 
synthetic-aperture  radars  (SARs)  by  means  of  convolution  with  the  Impulse 
response  of  a  two-dimensional  matched  filter,  is  given.  In  contrast  to 
other  theories  [1-14],  this  report  describes  a  two-dimensional  technique 
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that  applies  to  data  acquired  in  either  a  'squinted'  or  'sidelooking' 
configuration.  In  the  'sidelooking'  configuration  the  radar  antenna  is 
pointed  so  that  the  direction  of  the  centre  of  the  beam  is  perpendicular  to 
ground  track  [3],  In  the  'squinted'  configuration  the  antenna  points  in  a 
direction  other  than  the  perpendicular  [14]. 

In  order  to  convert  the  raw  radar  data  to  an  Image,  it  is  generally 
required  that  a  two-dimensional  signal  processing  operation  be  performed. 
Conventionally,  this  operation  is  segmented  into  two,  one-dimensional 
operations  where  the  radar  return  signals  that  are  associated  with  the  range 
(across-track)  and  azimuth  (along-track)  coordinates  are  independently  cross- 
correlated  with  their  respective  reference  functions.  Provided  the  range 
and  azimuth  signals  are  orthogonal,  and  provided  the  azimuth  signal  extent 
is  small,  this  approach  works  well. 

The  signal  processing  is  more  complicated  if  the  SAR  antenna  is 
'squinted',  or  if  the  azimuth  extent  of  the  signal  is  large.  In  these 
instances  the  range  and  azimuth  signals  are  coupled,  i.e.,  they  are  not 
independent  of  each  other,  and  the  operations  required  to  produce  a  high 
quality  distortion-free  image  are  more  complicated. 

In  the  following  sections  the  form  of  the  coupled  signal  is  examined, 
and  the  mathematical  operations  required  to  produce  a  high  quality  image  are 
described.  In  particular  a  novel  technique  that  employs  two-dimensional 
cross-correlation  of  the  received  SAR  signal  with  a  two-dimensional  reference 
function  is  described.  An  approximate  closed  formed  solution  is  given  for 
the  two-dimensional  correlation  integral  and  its  form  is  examined.  It  is 
shown  that  a  single  cross-correlation  function  can  be  used  to  produce  high 
quality  images  from  radar  signals  obtained  over  a  swath  in  range,  but  that 
these  images  must  be  geometrically  corrected  to  remove  positional  errors 
introduced  by  the  processing.  Equations  are  derived  to  describe  the  positional 
errors  and  the  operations  required  to  remove  them.  Finally,  the  modifications 
required  to  extend  the  theory  to  'multilook'  processing,  or  noncoherent 
averaging,  are  developed. 


2.  RANGE  AS  A  FUNCTION  OF  AZIMUTH  POSITION 

In  order  to  characterize  the  SAR  signal  that  results  when  the  trans¬ 
mitted  signal,  reflected  from  a  point  target  that  is  fixed  to  the  planet, 
is  returned  to  the  radar,  it  is  necessary  to  describe  the  range  to  the 
target  as  a  function  of  the  radar's  position  in  its  orbit.  To  do  this,  the 
flight  direction  of  the  radar  relative  to  the  surface  of  the  planet  is 
derived.  From  this  result,  an  equation  describing  the  range  to  the  point 
target  as  a  function  of  the  satellite's  orbital  position  is  obtained.  The 
satellite's  orbital  position  is  measured  relative  to  the  position  in  the 
orbit  at  which  the  point  target  is  in  the  centre  of  the  antenna's  horizontal 
pattern. 

The  velocities  of  the  vehicle  and  the  planet  are  calculated  In  the 
following  manner:  The  planet  is  assumed  to  be  a  sphere  of  radius  re  rotating 
with  angular  velocity  ue,  and  the  satellite  is  assumed  to  be  travelling  along 
a  circular  orbit  at  angular  velocity  us  (see  Figure  1).  The  tangential 


velocity  Vg,  of  the  reder,  with  respect  to  the  inertial  frame  is  u, (ra+hg) , 
where  hs  is  the  altitude  of  the  satellite  above  the  surface  of  the  planet. 

The  subsatellite  point  Si  (see  Figure  2),  is  the  point  of  Intersection  of  the 
surface  of  the  planet  ana  the  line  joining  the  radar  and  the  centre  of  the 
planet  (see  Figure  5).  The  point  Si  has  a  tangential  velocity  Vg/Ca,  with 


respect  to  the  inertial  fr 


where  the  constant  C.  is  defined  as 


C  -  1  +  -S 
a  r 


(1) 


The  surface  of  the  planet  is  moving  beneath  Si  at  a  tangential  velocity  V^, 
defined  with  respect  to  the  Inertial  frame 


Fl/unZ  Vthleh  Equhmtmt  VMocfty 
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Figure  3.  Radar  Geometry 


The  point  Sj  has  a  velocity  Ve_  with  respect  to  the  rotating  surface. 

_ _ 1  _  _  £  If  4^  — 


The  magnitude  of  Veq  is 


[v~+ft-%)2] 


and  its  direction  relative  to  V_/Ca  is  (see  Figure  2) 
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The  velocities  V  and  V  are  the  components  of  Ve  normal  to,  and  parallel 
to  Vg/Ca,  respectively.  epIn  Appendix  A,  Ven  and  Vgp  are  shown  to  be 

V  »  u  r  cos9.,  (4 

ep  e  e  i 


Ven  "  “ere  -in0iCO^o» 


where  6^  is  the  angle  of  inclination  of  the  orbit,  and  is  the  position  of 
the  spacecraft  in  the  orbital  plane  measured  from  the  equator  (see  Figure  1) 
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Equations  (2)  and  (3)  are  valid  for  the  vehicle  in  an  ascending  or  descending 
orbit. 


It  is  assumed  that  the  radar  antenna  is  pointed  at  an  angle  a  with 
respect  to  ground  track  as  shown  in  Figure  3.  In  this  report  the  angle  0y 
will  be  called  the  squint  angle.  The  angle 

(6) 

is  then  the  complement  of  the  squint  angle.  This  angle  is  composed  of  two 
components,  i.e.. 


6-0  +0 
y  ye  ya 


(7) 


where  0  is  caused  by  planet  rotation  and  6ya  is  the  antenna  pointing  angle 

relative  to  the  perpendicular  to  V  ,  The  angles  6ye  and  a  are  defined  with 

respect  to  the  equivalent  velocity  vector  Veq,  as  shown  in  Figure  2,  for  the 

case  6„  -0. 
ya 

The  slant  range  from  the  radar  located  at  S2»  to  a  point  target  A  on 
the  planet's  surface  (see  Figure  3),  is 

r2  *  Ire  +  (re+hs)2  -  2re^re+hs>CO802^JS*  (8) 


where  @2  is  the  angle  subtended  by  the  arc  AS2  on  the  surface  of  the  planet 
when  S2  is  the  sub-satellite  point. 


Equation  (8)  can  be  rewritten  to  show  how  r2  varies  with  both  the 
satellite's  orbital  position  and  the  complementary  squint  angle  a.  From 
spherical  trigonometry  and  Figure  3  it  can  be  shown  that 


CO80„  “  COS0  cos(6  +0.), 
2  r  a  b 


and 


sin6 

A 


cosasinO 
cos0f  * 


(9) 

(10) 


The  angles  8r,  0b  and  6  are  defined  as  follows)  6r  is  the  angle  subtended 
at  the  planet's  centre  by  the  arc  AS0,  where  S0  is  the  sub-satellite  point 
when  the  satellite  is  located  at  its  point  of  closest  approach  to  A;  6a  is 
the  angle  which  defines  the  satellite's  position  in  the  orbit  where  the 
radar  antenna  is  pointing  directly  at  the  target  at  A;  0^  is  the  angle 
measured  relative  to  0a  which  defines  the  present  position  of  the  satellite; 
and  6  is  the  angle  subtended  by  the  arc  ASj  on  the  surface  of  the  planet, 
when  S}  is  the  sub-satellite  point.  If  one  expands  cos(0a+6^)  in  (9)  and 
substitutes  for  sinOa,  by  using  (10),  the  equation  for  r«  can  be  rewritten 
as  a  function  of  0^. 

r2(0b)  -  |r2  +  (re+h#)2  -  2r#(re+ha>  [cos6eos6b  -  sin0coaasin0b]|,S. 
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It  is  more  convenient  to  work  with  r2(6b)  expressed  as  a  polynomial, 
than  with  the  form  given  by  (11).  The  polynomial  form  is  obtained  by 
expanding  rodv)  in  a  Mclaurln  series  about  6a.  Previous  treatments  [1]  have 
expanded  r2  (0^,;  about  0a*O.  An  expansion  about  0a“O  is  adequate  when  there 
is  no  squint,  i.e.,  a  ■  ir/2  rad;  however  for  a  a  r  r/2  rad.  an  expansion 
about  0a  as  given  by  (10),  yields  a  series  which  requires  fewer  terms  to 
represent  r2(0b)  accurately. 

A  Mclaurin  series  expansion  for  r2 *  about  angle  0a  is 


.  r  <p'(0) 

-  ;0  —  V- 

The  first  three  derivatives  of  the  series  evaluated  at  0^*0  are: 

r2 (°)  “  rl 

r  d)f(n  .  re(W  8iner 
r2  '  '  r,  tano 


r2<2)(0)  = 


2r, 


tl  +  (r+h) 2  r  [r  ^(0)] 

e _ e  S  1  l  _ 

r. 


and 


r2(3)(0)  -  -r2(1)(0) 


1  3  Tl  +  <r.+h.)2  3't2<1><0)|2 

‘  2  +  2  r2  ‘  r2 
'l  1 


(12) 

(13) 

(14) 

(15) 

(16) 


where  r^  is  the  slant  range  to  point  A  at  angle  0a.  It  can  be  shown  that  in 
most  cases  the  terms  in  0£  and  above  in  (12)  can  be  neglected.  Thus  r2 (0^) 
can  be  written 


r2(s)  “  aQ  +  ajS  +  a2s2, 

where  s  is  the  arc  length  along  the  ground  track  and  is  given  by 


and. 


8  “  ebV 


*o  “  V 


r  sine 
al  "  r^  Ca  tana  * 


(17) 

(18) 

(19) 


and. 


(20) 


f  a  ^ 

1  JL_  (1+cA  _  JL  „  JL_ 

2  2r.  U  a}  2r*  r, 

1  el 


Equation  (17)  is  the  general  expression  for  the  slant  range  to  a  point 
reflector. 

For  a=n/2  (17)  simplifies  to 

r2(s)  =  aQ  +  a2s2 

where 


2^2^  ^1+Ca^  "  2r|j  » 


and  rQ  is  the  slant  range  to  point  A  when  the  radar  is  at  its  point  of 
closest  approach  to  A. 

A  general  expression  for  range  as  a  function  of,  time  (or  orbital 
position  0^),  is  given  by  (17).  Equation  (17)  is  an  approximation  that  is 
valid  for  a  small  angular  extent  about  0a.  The  error  introduced  by  approxi¬ 
mating  (11)  by  (17)  can  be  evaluated  by  calculating  the  difference  between 
these  two  equations  for  the  values  of  interest. 


In  the  derivation  of  (17)  a  circular  orbit  and  a  spherical  planet, 
i.e.,  constant  hs  were  assumed.  This  constraint  can  be  relaxed  so  that  an 
elliptical  orbit  and  ellipsoidal  planet  can  be  accommodated .  Appendix  B 
gives  the  derivation  of  an  extended  version  of  (11)  which  includes  a  constant 
vertical  velocity  Vv.  Obviously,  for  the  more  general  case  of  an  elliptical 
orbit,  neither  Ve<j  norVv  is  truly  constant.  However  for  realistic  orbits, 
accelerations  $  and  Vv  are  small. 


3.  THE  FORM  OF  THE  TWO-DIMENSIONAL  SIGNAL  FOR  A  POINT  TARGET 

In  this  treatment  of  SAR  it  is  assumed  that  the  surface  causing  the 
radar  backscatter  can  be  modelled  by  a  collection  of  point  scatterers.  The 
SAR  illuminates  and  receives  the  signal  scattered  from  these  scatterers.  In 
this  section  this  signal  is  described.  Subsequent  sections  will  discuss  the 
signal  processing  that  operates  upon  the  received  signal  from  each  point 
scatterer  and  produces  an  image  of  each  point. 

The  term  twc-dimenslonal  is  used  to  describe  the  signal  because  the 
received  signal  is  formed  by  gathering  information  from  each  pulse  trans¬ 
mitted  (the  range  distension),  and  from  changes  that  occur  in  the  received 


signal  from  pulse  to  pulse  (the  along  track  or  azimuth  dimension),  as  a 
result  of  the  relative  motion  between  the  radar  and  the  scatterer.  The 
derivation  of  the  form  of  the  two-dimensional  signal  consists  of  the  follow¬ 
ing  steps: 

1.  the  determination  of  the  form  of  the  transmitted  signal; 

2.  the  derivation  of  the  form  of  the  signal  received  by  the  radar  when 
this  transmitted  signal  is  reflected  back  to  the  radar  by  a  point 
target ;  and 

3.  the  derivation  of  the  form  of  the  signal  obtained  by  mixing  in  a 
quadrature  detector,  the  received  signal  with  the  carrier  of  the 
transmitted  signal. 

It  will  be  described  in  later  sections  how  this  two-dimensional  signal 
is  processed  to  form  the  two-dimensional  radar  image  of  the  point  target • 


3.1  FORM  OF  THE  TRANSMITTED  SIGNAL 

As  the  vehicle  upon  which  the  radar  is  mounted  proceeds  along  its 
flight  path,  the  radar  transmits  a  series  of  pulses  (see  Figure  4)  of  the 
form 


fT(t> 


°°  j 

Rft  E  W(t)^(t-mT)exptj(«ct+<|>o)n  , 

ma— *  ) 


(26) 


where  \^(t-mT)  is  a  complex  modulation  function,  W(t)  is  a  weighting  function 
which  will  be  defined  later,  t  is  time,  T  is  the  interpulse  period,  m  is  the 
pulse  number,  $0  is  the  phase  of  the  carrier  at  t«0,  and  wc  is  the  carrier 
frequency.  The  argument  (t-mT)  of  the  modulation  function  represents  time 
measured  from  the  start  of  the  pulse.  It  is  assumed  for  convenience 
that  m*0  when  the  satellite  is  located  at  (see  Figure  3),  i.e.,  at  an 
angle  6a  from  the  position  occupied  by  the  satellite  when  it  is  at  the  point 
of  closest  approach  Sq. 


Figure  4.  Trenmitfed  Signet 
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3.2  FORM  OF  THE  RECEIVED  SIGNAL 

The  signal  transmitted  at  time  t  from  position  S^.  (Figure  11) ,  and 
received  at  time  t+At  at  position  SR,  after  reflection  from  the  point  target 
at  A,  is  given  by 

fR(t)  *■  Re |rw2 ( t ) i|»( t-mT- A t ) 

exp|j[o)c(t-At)  +  (27) 

where.  At  is  the  time  taken  by  the  signal  to  travel  the  path  S<{«ASr  (Figure 
11),  and  r  is  the  reflection  coefficient  of  the  point  scatterer  at  A.  It  is 
assumed  that  the  point  target  at  A  is  aspect  insensitive  so  that  T  does  not 
change  during  the  time  the  radar  views  the  target. 

The  range  to  the  target  changes  during  the  time  taken  to  travel  the 
path  SjASr  (Figure  11).  It  is  shown  in  Appendix  C  that  a  very  good  approxi¬ 
mation  for  At  is 

2r. (s) 

At  =  -  .  (28) 

c 


3.3  FORM  OF  THE  SIGNAL  OUTPUT  BY  THE  QUADRATURE  DETECTOR 

The  received  signal  fR(t)  is  demodulated  in  a  quadrature  detector  by 
mixing  fR(t)  with  the  carrier  of  fT(t+At).  The  transit  time  At  can  be 
replaced  by  2(a0+a1s+a2s2)/c  using  (17)  and  (28),  and  then  the  complex  signal 
output  by  the  quadrature  detector  can  be  expressed  as  the  two-dimensional 
function 

fD(t’,s)  -  j  W2(t,,s)<|»^t'  (aQ4«18+a2s2)J  exp[-j2k(a£)+a1s+a2s2)],  (29) 

where. 


t*  «  t-mT,  (30) 

is  time  measured  from  the  start  of  the  pulse, 

s  -  mVe<T,  (31) 

is  the  arc  length  travelled  by  the  sub-satellite  point  along  the  sub-satellite 
track  during  m  interpulse  periods,  and 

k  -  ^  ,  (32) 

is  the  wave  number.  The  variable  t'  is  a  measure  of  slant  range.  The 
specific  value  of  t'  given  by  At  (28)  refers  to  the  transmit  time  for  a 
target  at  a  specific  slant  range  r2> 
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In  eqn,  (29)  W2(t)  has  been  rewritten  as  the  two-dimensional  function 
W(t*,s)  by  segmenting  W(t)  into  consecutive  sections,  of  length  T,  and  using 
the  variable  vsv  to  index  the  sections. 

The  weighting  function  W2(t',s)  does  two  things:  it  limits  the  range 
extent  of  the  strip  of  terrain  from  which  signals  are  received;  and  it  limits 
the  number  of  pulses  which  can  be  coherently  integrated  to  form  the  SAR  image 
of  the  target  at  A,  In  limiting  the  number  of  pulses  that  can  be  coherently 
integrated,  W2(t*,s)  also  limits  the  maximum  azimuth  resolution  attainable. 

It  is  assumed  that  the  position  and  width  of  W2(tf,s)  in  the  t1 
dimension  are  chosen  so  that  the  radar  illuminates  only  a  narrow  strip  of 
the  planet *s  surface,  containing  targets  at  slant  ranges  in  the  interval 

|  nT  <  r2  <  j  (trt-l)T,  (33) 

where  n  is  an  integer.  There  is  a  range  ambiguity,  therefore  the  value  of  n 
must  be  known  to  convert  At  correctly  into  r2  via  (28) • 

When  W2(t(,s)  is  chosen  to  satisfy  (33),  signals  at  ranges 

r’  =  r2  ±  P  f  T,  P  -  1,2 .  (34) 

which  would  otherwise  arrive  at  the  radar  during  the  time  interval 
(nT  t*  (rH-l)T]  and  produce  ambiguous  signals  superimposed  on  the  desired 
signal,  do  not  have  to  be  considered.  It  is  also  assumed  that  the  Interval 
of  r 2  which  satisfies  (33)  is  large  compared  to  the  range  extent  in  t*  of  the 
modulation  function  ip  in  (29),  and  that  W2(t,,s)  varies  slowly  with  tf.  With 
these  assumptions  W2(tf,s)  can  be  written  as 

W2(t',8)  -  WqW^s)  ,  (35) 

where  is  a  constant. 

The  shape  of  U^(s)  in  the  s- dimension  is  assumed  for  ease  of  analysis 
to  be  rectangular.  The  actual  shape  of  W^(s)  is  determined  by  the  antenna 
pattern.  The  following  analysis  will  differ  in  the  details  but  not  the 
general  concept  if  a  different  W^(s)  is  used.  Equation  (34)  for  W2(t',s), 
can  be  written 

W2(t',s)  =  WQ  rect^-^  ,  (36) 

where  Lg  is  the  synthetic-aperture  length.  The  size  of  Lg  determines  the 
number  of  pulses ,  returned  from  the  point  target  at  A,  which  can  be  coherently 
processed  to  form  the  SAR  image  of  A. 

The  general  two-dimensional  signal  for  a  unity  gain  point  target  is 
obtained  by  setting  r  to  unity  in  (29),  i.e., 

h1(t* ,8)  -  WQ  rect/~-^<J»|t’  -  ~  <a0+al8+a282>]  exp [ - j 2k (a^a^a^ 2 ) ] .  (37) 
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Figure  5  shows  the  form  of  h^(t'(s).  The  responses  from  three  point 
targets  are  shown.  At  some  particular  instant  in  time  targets  (1)  and  (2) 
both  lie  in  the  same  direction  a  but  are  at  different  rj/s.  Target  (3)  lies 
at  the  same  range  rg,  as  target  (1),  but  is  situated  at  a  different  o. 

The  thick  solid  lines  in  Figure  5  represent  the  envelope  of  h^(t',s). 
The  lines  limiting  the  s  extent  are  the  edges  of  the  rect(s/Lg)  function. 
This  function  increases  in  length  as  the  radar's  antenna  pattern  broadens 
with  range.  In  other  words  the  arc  length  Lg  over  which  the  radar  views  the 
target  increases  with  range. 

The  curved  lines  limiting  the  t'  extent  of  the  signal  are  the  edges  of 
|i|/|  in  (37).  The  curvature  of  these  lines  is  determined  by  the  manner  in 
which  r2  varies  with  s  (shown  by  the  dashed  lines  in  Figure  5). 

The  thin  solid  lines  inside  the  ({/-envelope  represent  constant  phase 
contours  of  <{/.  The  phase  of  t{/  is  composed  of  two  factors:  the  phase 
exp[~j2k(ag+ais+a2S2) ]  along  the  's'  coordinate  (see  Figure  5),  and  the 
phase  of  <J/  along  the  t*  coordinate. 


4.  SIGNAL  COMPRESSION 

The  signal  (29)  returned  from  the  point  target  at  A  (Figure  3)  must  be 
processed  to  produce  a  two-dimensional  image  of  the  target.  The  goal  of  this 
section  is  to  show  the  form  of  the  solution  resulting  from  the  two-dimensional 
convolution  of  the  signal  received  from  a  point  target  (r*l)  with  the  two- 
dimensional  impulse  response  of  a  filter  that  is  matched  to  the  signal  at  a 
particular  reference  range.  A  mathematical  analysis,  that  results  in  an 
approximate  closed-form  solution  for  the  equation  describing  the  image,  is 
presented. 


Figure  S.  The  form  of  the  tignalt  received  from  three  point  target*.  Target*  (1)  and  (21  are  Illuminated  by 
the  antenna  at  the  tame  time.  Target  131  it  at  the  tame  rQ  at  target  (!)  but  I*  illuminated  at  a  later  time. 


In  performing  the  analysis  the  general  point  target  is  assumed  to  be  at 
range  r1  and  complementary  squint  angle  of  and  to  have  an  associated  synthetic- 
aperture  length  Lg.  the  two-dimensional  Impulse  response  of  the  matched  filter 
is  calculated  for  a  point  target  with  a  reflection  coefficient  of  unity, 
located  at  a  reference  range  r^,_and  complementary  squint  angle  a,  and  which 
has  a  synthetic  aperture  length  Ls.  _ The _ values  rj,  a  and  Lg  are  chosen  to  be 
different  from  the  reference  values  r^,  a  and  Lg.  This  is  done  so  that  the 
effect  of  mismatch  between  the  signal  and  reference  function  can  be  examined 
in  the  final  closed  form  solution  for  the  SAR  image. 

4.1  THE  TWO-DIMENSIONAL  MATCHED  FILTER 

The  impulse  response  of  the  matched  filter  is  the  time-reversed  conju¬ 
gate  of  the  two-dimensional  signal  returned  from  a  point  target  (37).  When 
the  parameters  used  in  (37)  are  changed  to  the  reference  values,  the  impulse 
response  of  the  matched  filter  can  be  written  as 


h2(t',8)  -  h*(-t’,-s) 


l.e.. 


h2(t',s)  -  recti 


"  «  (Vai8+i2#2)]  •*PtJ2k<Vai^i2,2)1»  (39) 


where  the  a q,  3^,  and  &2  refer  to  the  reference  values  of  sq,  a^,  and 
which  are  found  by  substituting  and  5  in  (19)-(21). 

In  the  actual  analysis  a  slightly  modified  form  of  (39),  in  which  the 
3q  term  has  been  removed,  will  be  used 

h^t'.s)  -  rect^j^ ^  (-ais+a2s2)j  exp[j2k(-l1»t-a2s2)].  (40) 

The  use  of  hoCt' ,s)  instead  of  h2(t',s)  causes  targets  to  be  mapped  at  range 
aQ  instead  of  range  a q-Zq  in  the  final  SAR  image.  In  other  words  a  shift  of 
'the  map  origin  is  introduced. 

4.2  SOLUTION  OF  THE  TWO-DIMENSIONAL  CONVOLUTION  INTEGRAL 

An  equation  describing  the  form  of  the  output  SAR  image  is  obtained 
by  solving  the  two-dimensional  convolution  Integral 

°o  wr-  /  .  .  v 

f0(t',s)  -JJ  h3(t'-n,s-C)fD(n,C)dOdC  r*Ct\^/r*Ct\L~) 

exp|-j2kla0+a1t-ai(?-s)  + 

/  “  I*1(C-»)  ♦  a2(C-.)2]}  p[q  -  |  (a^frfa^2)!  dfldC.  (41) 
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The  solution  to  (41)  will  he  obtained  in  three  steps.  The  first  step 
will  be  to  solve  the  Inside  integral  by  Integrating  over  fl.  The  second  step 
will  be  to  place  a  restriction  on  the  £  variation  of  the  solution  to  the 
inside  integral.  Provided  certain  conditions  are  satisfied  this  restriction 
will  allow  a  separation  of  the  variables.  The  third  step  will  be  to  solve 
the  outside  integral  which  is  a  function  of  (. 

4.2.1  Solution  of  the  Inside  Integral 

The  solution  to  the  inside  integral  is  g(y)  where 

Y  “  t'  -  ~  [a^a^-S^C-aHa^-g^C-a)2]  (42) 

The  function  g(y)  is  the  auto-correlation  of  p,  and  has  its  maximum  at 

t’  “  f  [a0+a1C-ai(C-s)+a2C2-*2a-s)2l.  (43) 

The  form  of  the  modulation  function  ^  has  not  yet  been  specified. 
Modulation  functions  used  for  pulse-compression  normally  have  the  character¬ 
istic  that  most  of  the  energy  in  the  output  pulse  is  confined  to  a  main  peak 
which  has  a  3  dB  width  approximately  equal  to  irc/Kx,  where  2ir/Kr  is  the  band¬ 
width  in  Hz.  In  this  report  only  modulation  functions  of  this  type  will  be 
considered. 

The  shape  of  the  g(y)  function  that  is  obtained  for  a  typical  modulation 
code  <J»(t')  is  examined  next.  If  a  linear  FM  modulation  function  is  assumed, 
4*(t')  is  given  by 

<Kt')  =  rect^ - exp  |j  |  ^t*  -  y ^  , 

where  t  is  the  pulsewldth,  and  K  is  the  linear  FM  rate  in  rad/s2, 
given  by  (44) , 

g(y)  ■  (*-|y|)  rect^f^  sine  |^y  (x-|y|)yJ  . 

where  sine  x  =  sin  x/x. 

The  shape  of  g(y)  is  shown  in  Figure  6.  The  sinc(KY/2  (t— | y | ) 1  factor 
produces  three  peaks;  the  outer  two  of  which  are  suppressed  by  the  (t-|y|) 
factor.  Provided  the  time-bandwidth  product  of  the  range  (t*)  signal  is 
high,  i.e.,  provided  y  «  t  In  the  vicinity  of  Y“0»  the  major  portion  of  the 
signal  energy  will  be  contained  in  the  central  peak.  The  equation  describing 
the  form  of  this  peak  can  be  closely  approximated  by 

g(Y)  “  t  sinewy  ty^  .  (46) 

The  r,  resolution  is  proportional  to  the  malnlobe  width  (the  width  between 
the  first  null  on  each  side  of  y*0),  which  is  given  by 


(44) 

For  p(t') 

(45) 
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Figure  6.  Form  of  Signal  from  a  Single  Point  Target.  After  the  Solution  of  the  Inside  Integral. 


R 


2itc 

Kt 


(47) 


It  can  be  seen  in  Figure  6  that  the  position  of  the  central  peak  produced  by 
the  solution  of  the  inner  integral  is  at  Y“0.  Froa  (42)  It  is  known  that  y 
depends  on  ^"s)"a2(Z_*2 21  which  in  turn,  is  a  function  of  the 

azimuth  (£)  position.  The  C -dependence  is  a  manifestation  of  the  coupling 
or  interdependence  of  the  range  (t’)  and  along-track  (a)  signals. 


4.2.2  Separation  of  Variables 


In  order  to  obtain  an  analytical  solution  to  the  outside  integral  in 
(41).  y  in  g(y)  must  be  constrained  such  that  the  ^-dependence  is  removed. 
If  one  lets 


A 


2 


a2“82 


(48) 


arai+2V 


and 


(49) 


V41 


2j2* 


(50) 


one  can  rewrite  (42)  as 

y  -  f  -f  (AjC2^!^).  (51) 

The  £ -dependence  can  be  removed  from  (46)  by  choosing  A2  and  At  in  such  a  way 
that  the  position  of  the  main  peak  does  not  deviate  from  y“0,  by  more  than  a 
■•11  fraction  of  its  width;  l.e. ,  a  small  fraction  of  2vc/Kt.  This  deviation 
is  described  by  the  following  equation 
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«r  '  'V2  +  V'  «  2  ft* 

With  the  restriction  given  by  (52),  the  substitution 

Y  »  t«  -  \  (Aq  ±  «r) 

can  be  made  for  y  in  g(Y>.  and  an  approximate  analytical  solution  for  (41) 
obtained. 

4.2.3  Solution  of  the  Outside  Integral 

The  remaining  step  in  the  solution  of  (41)  is  the  integration  over  5. 
With  y  given  by  (53),  (41)  can  be  written  as 

C-H* 

fp<f,s)  -  -f- *['■  -  c  <W]  /  r*Ct(^i) 

rect^-^  exp[-j2k(Ag  +  A^£  +  A^S2) ]d£  •  (54) 

The  solution  to  (54)  has  two  major  cases,  A2*0  and  A2/O,  and  two  sub¬ 
cases  of  each  of  the  major  ones,  Lg  >_  Lg  and  Lg  <  Lg. 

The  variable  A2  has  a  special  significance  in  that  it  determines  the 
magnitude  of  the  quadratic  phase  error,  i.e., 

*e  -  IkA^2.  (55) 

As  shown  in  [15]  the  magnitude  of  the  quadratic  phase  error  determines  whether 
or  not  an  image  can  be  considered  focussed.  In  [15]  the  in-focus  condition  is 
defined  by 


(52) 

(53) 


|*J<J*md  .  (56) 

In  this  section  the  solution  for  A2"0  will  be  obtained,  and  then  by 
reference  to  [15]  it  will  be  shown  that  for  the  solution  gradually 

degrades  as  | A2 1  Increases.  With  this  link  established,  it  will  be  assumed 
that  the  solution  for  A2m0  can  be  used  as  long  as  the  image  is  in  focus. 

For  Lg>Lg,  the  impulse  response  of  the  matched  filter  is  longer  than 
the  signal  from  which  the  synthetic-aperture  is  formed.  It  will  be  shown 
later  that  in  certain  situations  Lg  must  be  greater  than  Lg  in  order  to 
obtain  the  full  along-track  resolution. 

For  LS<LS,  the  impulse  response  of  the  matched  filter  la  shorter  than 
the  signal  from  which  the  synthetic-aperture  is  formed.  Sometimes  Lg  is 
chosen  to  be  smaller  than  Lg  because  the  antenna  along-track  illumination  is 
broader,  and  thus  Lg  is  longer,  than  required  to  obtain  the  desired  resolu¬ 
tion. 
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4. 2. 3.1  Solution  for  Zero  Quadratic  Phase  Error 

u  ls  -  Ls 

For  this  sub-case  the  solution  to  (54)  is  given  by 

fo(t',S)  =■  g|t'  -  |  (A0±6r)Jexp(-j2kA0) 


/  L  \ 

vl 

r  +  2)/  W\ 

jkAj^ 

-  -  Ml 

recti  1  Is  +  — 2 — 1  exP 

*  ! 

\  s  /  ' 

L  J 

i 

sine 


+r,ct[w] 


L  J 

f(.4)l 

Lg  sine (kA^Lg)  +  rect  L 

_  S  • 


(-8  +  -  V1)  exp  jkAi  ("8  +  "V2) 


8  Inc 


(57) 


Figure  7  shows  the  form  of  (57).  It  can  be  seen  that  it  is  highly  desirable 
to  ensure  that  the  position  of  the  main  peak  of  the  azimuth  response  Sp, 
where 


Vai 


p  2k 


is  confined  to  the  interval 


<VLS> 


2 

(L.-LJ 


<  8  < 


2  -  p  -  2 


(58) 


(59) 


In  this  Interval  neither  attenuation  nor  broadening  of  the  main  peak  occurs, 
and  the  main  peak  width  (i.e.,  between  the  nulls),  la  given  by 


(60) 


Figure  7.  Form  of  Signal  From  a  Single  Point  Target,  After  the  Solution  of  the  Outside  Integral 
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For  (59)  to  hold,  Lg  must  be  chosen  such  that 


s-  Ls  ‘  Ls 


Ls  -  Ls  +  “p‘ 


For  this  sub-case  the  solution  to  (54)  Is  given  by 

f  (t’,s)  *|(vv] 


(62) 

(63) 


the  naln  peak  width,  between  the  nulls,  ls 


4. 2. 3. 2  Solution  for  Non-Zero  Quadratic  Phase  Error 


The  purpose  of  this  section  is  to  obtain  the  solution  for  (54),  when 
A2  J*  0.  By  reference  to  [15]  it  is  shown  that  as  A2  approaches  zero,  the 
solution  to  (54)  gradually  assumes  the  form  of  (57)  or  (62). 

Uhen  A2  0,  the  solution  to  (54)  is  found  by  solving  the  following 
integral 


The  solution  to  (65)  has  two  sub-cases,  Lg>L§  and  Lg<Lg.  A_link  with  (57) 
and  (62)  will  be  established  by  examining  the  solution  for  Lg>Ls.  The  solu¬ 
tion  for  Lg<Lg  is  found  in  a  similar  manner.  When  Lg>Lg,  the  solution  to 
(65)  is 


+  rect 


jtc<«4>  -  c«3>  +  ^s«4>  -  S«3>1 


+  rect 


[c(5l)  -  C(5')]  +  J[S«»)  -  8(5 


;>i||  . 


(61 


20 


are  Freanel  Integrals,  and 


The  second  term  In  (66)  represents  the  solution  when  the  main  peak  is 
located  in  the  region  (Lg-Lg)/2  to  (Lg-Lg)/2.  In  this  region  (66)  becomes 

r  (f..)  -  ^  g(r)  «p  -32k  (*„  - 

2  (nj  {‘C<P>  *  C«»  +  J'S<P>  *  S«»}  •  <73 


where,  using  the  notation  of  [15], 


« -7=  (£-“)-<;• 


I,  -  Ls  ^  -  /5T 


The  substitution  in  terms  of  At*  and  n  allows  comparison  of  (73)  with  the 
results  in  [15].  In  [15],  it  is  shown  that  as  ♦smf»  decreases  from  w/8  rad 
to  0  rad  only  a  slight  change  in  the  shape  of  the  main  compression  peak 
occurs.  Similarly,  the  solutions  (57)  and  (62)  degrade  only  slightly  as 
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| A2 1  ■*  n/8  rad.  Thus,  for  ♦ewav,  <_  ir/8  rad  (57)  and  (62)  can  be  used  as  the 
general  solution  to  (54). 

4.3  PROCESSING  SWATH  IN  RANGE 

The  allowable  processing  swath  in  range  Ar^  over  which  a  single  matched 
filter  can  be  used,  it  determined  by  the  maximum  range  mismatch  which  satis¬ 
fies  both  the  restriction  on  $e,  (56),  and  the  equation  for  6r,  (52). 

4.3.1  Restriction 

As  can  be  seen  from  (55),  the  maximum  quadratic  phase  error  occurs  at 

£  which  from  (54)  is 
max 


The  maximum  value  of 


Cm ax  “  2 


Arl  *  2 1 rl"f 1 1  » 


is  found  by  solution  of  the  following  equation: 


(r)  if  • 


where  4>e  has  been  replaced  by  (55),  and  £  by  (78).  After  the  substitution  of 
(48)  for  A2,  and  of  (32)  for  k,  (80)  can  be  written  as 


8L2 

—  Ia2“a2^  —  1  * 


After  replacement  of  a2  and  I2  in  (81),  the  equation  can  be  rewritten  as 
8L2  I  -  fl  +  C2  r2  r2  sin20  “I 

_S  JL  _ a  _  _i _ •  C2  _£  I 

X  2^  ^  2  2r2  r£  a  tan2a  J 

-  [*1  +  C2  f2  r2  sln20  *11 
"  2f^  [  2  *  “  2r2  “  ff  Cl  tan^JI  £  1.  (8J 

If  it  is  assumed  that  (r2/r2)  «  1,  and  a-8,  then  (82)  can  be  simplified  to 


BSL, ; 
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4.3.2  6r  Restriction 

In  this  section  it  is  shown  that  a  suitable  6r  can  be  found  to  satisfy 
(52).  Bounds  have  already  been  placed  on  £  (78),  and  A2  (80).  Thus  the  only 
step  remaining  is  the  derivation  of  a  bound  on  A^.  The  validity  of  the 
solution  to  (40)  is  restricted  to  a  neighborhood  of  radius  2  Waz  about  the 
main  compresion  peak.  With  this  restriction,  |AjJ  is  limited  to  the  Interval 

lAll  <  K  -  ai  +  2a2(8p  *  2  WaZ) I ‘  (88) 

After  substitution  for  sD  (59)*  and  W  ((60)  or  (64))*  (88)  can  be  written 

r  az 

IaJiN.  (89) 

If  the  restriction  (89)  is  coupled  with  those  on  A2  (80)  and  £  (78) *  then 
(52)  can  be  rewritten  as 


or 

6  = 

r 

Equations  (90)  and  (91)  give  the  maximum  value  of  6r*  The  width  Wr  is  always 
much  greater  than  A  because  of  fundamental  physical  constraints.  Therefore 
the  solutions  given  by  (57)  and  (62)  are  only  constrained  by  the  quadratic 
phase  error  restriction  (56).  Henceforth*  the  quadratic  phase  error  restric¬ 
tion  will  be  considered  as  the  only  restriction. 

The  restriction  on  5r  ensures  that  no  additional  energy  is  spilled 
into  the  two-dimensional  sidelobe  structure  by  the  range-azimuth  coupling. 

The  solutions  obtained  subject  to  this  restriction  do  not  provide  any  insight 
as  to  the  detailed  form  of  this  sidelobe  structure.  The  sidelobes  can  only 
be  fully  described  by  a  numerical  solution  of  (41). 


—  ±  \ 
32  A 


« 


V  Ls  -  Ls» 


(90) 


xli 
32  Ls 


«  Wr,  Ls  <  Ls 


(91) 


4.4  PROCESSING  INTERVAL  IN  AZIMUTH 

In  the  along-track  or  azimuth  dimension,  there  are  also  restrictions 
on  the  processing  swath.  As  the  vehicle  carrying  the  radar  moves  along  its 
track  it  will,  in  general,  traverse  lines  of  constant  latitude.  The  latitude 
change  introduces  a  change  in  0ye,  and  therefore  a.  If  a  changes  appreciably 
in  value  from  5,  the  reference  function  (40)  will  no  longer  match  the 
properties  of  the  signal  received  by  the  antenna.  The  length  of  the  refer¬ 
ence  function  Lg  must  be  chosen  according  to  (63)  so  that  it  ls  long  enough 
to  properly  match  the  data  that  ls  received. 

4.5  POSITIONAL  ERRORS  IN  COMPRESSED  SIGNAL 

In  this  section  It  is  shown  that  the  two-dimensional  convolution 
selects  portions  of  the  signal  fp  with  slope  S^.  These  selected  portions  - 
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are  shown  to  be  offset  In  and  s  from  the  centre  of  the  two-dimensional 
signal  returned  from  the  target. 

The  selection  process  can  be  understood  by  examining  Figure  8.  In 
Figure  8  are  shown  three  sets  of  lines:  the  dotted  is  a  line  of  constant  gjj 
the  solid  are  lines  of  constant  S|  and  the  dashed  are  lines  of  constant  rQ. 
The*  slopes  of  the  dashed  lines  are  given  by 


dT  “  *1  +  2a28  * 

For  the  particular  line  through  D,  (92)  can  be  written  as 


3 r  -  *i-  <”> 

Equation  (93)  is  the  slope  of  the  phase  fronts  of  the  two-dimensional  refer¬ 
ence  function  ho (t ' ,o) .  The  two-dimensional  convolution  will  cause  a  section 
of  fp(t',s)  to  be  selected  that  has  a  slope  g^.  It  is  known  from  (64)  that 
the  compression  peak,  for  a  target  located  along  the  line  s-0  (e.g.,  E  in 
Figure  8),  will  be  located  at  s_  (e.g.,  F  in  Figure  8).  For  &2  ~  *2  t^e 
slope  at  F  is 


where  df2/d*|  „q  i*  the  slope  at  D,  the  point  to  which  the  reference  function 
is  matched.  8Sxnce  the  locations  of  the  points  E  and  F  are  general,  it  can 
be  stated  that  the  convolution  does  indeed  select  portions  of  fD(t* ,s)  with 
the  same  slope  as  h^t'ys)  at  D. 


1st 


I 
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Ffgun  8.  Form  of  the  Itodopr 
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In  Figure  8  it  Is  shown  that  the  output  signal  is  also  shifted  in  range 
by  the  displacement  from  E  to  F.  The  output  position,  as  seen  in  (53),  is 
shifted  from  to  AQ  t  6ft  where 

A0  "  *0  +  Ar(V»  (95) 

and 

Ar(s  )  -is-  I,a*.  (96) 

P  1  P  2  P 

Both  the  r^-shift  Ar (s_) ,  and  the  s-shlft  sp,  are  range  dependent* 

They  must  be  removed  in  order  to  form  a  geometrically  correct  image. 


5.  IMAGE  FORMATION 

In  order  to  form  a  geometrically  correct  image,  the  azimuth  offsets  s 
(59),  and  the  range  offsets  Ar(s_)  (96),  must  be  removed.  After  removal  ofp 
these  offsets,  a  line  of  points,  e.g.,  AS£  in  Figure  9,  will  lie  along  a  line 
of  constant  s  in  the  (r^,s)  coordinate  system,  l.e.,  will  be  plotted  perpen¬ 
dicular  to  s.  When  viewed  from  the  (d0,s0)  coordinate  system,  where  (d0,s0) 
forms  an  orthongonsl  map  grid  on  the  surface  of  the  planet,  the  axes  and 
s  are  neither  coplsnar  nor  orthogonal.  Therefore  when  an  orthogonal  map  is 
required,  images  produced  in  the  (r^,a)  coordinate  system  must  undergo  a 
coordinate  transformation.  This  transformation  Involves  two  operations: 

1.  conversion  from  slant  range  r^  to  ground  range  d^,  and 

2.  conversion  from  processor  coordinates  (d^,a)  to  map  coordinates 
(do"*o)- 

Ground  range  is  obtained  from  slant  range  by  using  the  equation 


The  conversion  from  (d^,s)  to  (d0,So)  is  made  by  writing  the  constant 
s  line  (AS*  in  Figure  9)  at  an  angle  with  reapect  to  the  sub-satellite  track 
(80si).  The  data  must  be  laid  down  on  the  map  along  the  llnea  described  by 
the  intersection  of  the  antenna  centre-beam  plane  and  the  surface  of  the 
earth,  when  the  satellite  is  at  orbital  position  6a(ra+hs)«  The  equations 
which  describe  the  transformation  are 


and 


(98) 

(99) 
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6.  MULTI-LOOK  SIGNAL  PROCESSING 

This  section  describes  aulti-look  signal  processing*  e  technique  that  is 
employed  in  SAR  to  reduce  target  scintillation. 

In  Section  3.3  it  was  shown,  (see  eqn.  64))  that  the  azimuthal  resolu¬ 
tion  is  proportional  to  the  length  of  the  signal  function  that  is  used  to 
form  the  synthetic-aperture.  Frequently,  the  antenna  beamwidth  B  is  suffi¬ 
ciently  broad  that  the  resulting  signal  function  is  longer  than  that  required 
to  produce  the  desired  azimuth  resolution.  Rather  than  truncating  the  signal 
and  thereby  wasting  the  energy  associated  with  this  additional  signal,  it  is 
common  practice  to  segment  it  into  several  pieces,  each  of  appropriate  length 
for  the  resolution  desired.  Each  of  these  segments  is  matched- filtered  to 
produce  an  image. 

The  images  are  geometrically  corrected  and  converted  to  map  coordinates. 
After  registering  the  images  with  respect  to  one  another,  they  are  noncoher- 
ently  summed  together  to  form  a  composite  image.  The  terms  multi-looking, 
angle  diversity,  and  mixed  integration  are  the  jargon  commonly  used  to 
describe  the  noncoherent  averaging  operation. 

The  segmenting  operation  described  above  results  in  using  portions  of 
the  azimuth  signal  that  are  acquired  at  different  o-angles ,  (see  Figure  9). 
Each  of  these  images  are  referred  to  as  "looks";  the  qth  look  being  at  angle 
a  .  Obviously  the  a. 's  must  be  selected  such  that  the  corresponding  signal 
segments  are  contained  within  the  antenna  horizontal  beamwidth. 


In  order  to  register  the  looks,  it  is  necessary  to  specify  the  matched 
filters  for  each  of  the  looks  so  that  each  filter  corresponds  to  the  same 
point  A.  To  specify  the  qth  filter,  a  la  first  selected,  subject  to  the 
constraint  mentioned  above,  and  Sq.,  a]L,  and  Sj-  are  calculated  using  riq 
and  ciq.  The  parameter  r*q  is  the  slant  range  toHA  when  A  is  at  an  angle-to- 
track  oq.  In  terms  of  the  fundamental  parameters  a  and  rj,  r^q  is  given  by 


The  q 


th 


filter  is  then  specified  by 


*0q 


(101) 


(102) 
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(103) 


(104) 


The  matched  filtering  and  image  formation  operatlona  for  the  multi-look 
case,  are  Identical  to  those  required  for  the  single  look  case,  except  that 
In  all  the  equations  the  parameters  a.,  a.  ,  a.  ,  a  end  r.  are  substituted 
for  aQ,  a^  a2,  o,  and  r^.  q  lq  q  q  Aq 


7.  EXAMPLE 

A  software  digital  processor  based  on  the  theory  described  In  this 
report  has  been  programed  and  used  to  produce  Imagery  from  data  acquired 
with  airborne  and  satellite-borne  SARs.  Examples  of  processed  Imagery  are 
shown  In  Figures  12  and  13*  The  Images  were  produced  using  data  from  the 
SEASAT-A  satellite.  The  image  in  Figure  12  Is  a  four-look  25m  ground  range  x 
25m  azimuth  resolution  image  of  Trols  Rivieres (  Quebec.  The  Image  In  Figure 
13  is  a  single-look  25m  ground  range  x  7m  azimuth  resolution  image  of 
Halifax,  Nova  Scotia.  The  processing  parameters  for  these  two  scenes  are 
given  in  Table  1.  For  the  SEASAT-A  satellite  the  antenna  beam  is  only 
slightly  squinted  (3-4*);  however,  the  relatively  low  frequency  (1.276  GHz) 
of  the  transmitted  signal,  and  the  large  distance  between  the  satellite  and 
the  planet's  surface,  mean  that  procesalng  the  data  to  7  and  25m  resolutions 
involves  the  same  degree  of  complexity  as  more  highly  squinted  airborne 
systems. 


Fl&nfl.  DtfkHtkm  of  Slant  tongmst  Ttwrnnftmd  NsesA*  Hum 


4  LOOKS 


Figure  12.  SEASA  T-A  SAR  Image  of  Troit  Rivieres,  Quebec 
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Figure  13.  SEASA  T-A  SAR  Image  of  Halifax ,  Nova  Scotia 
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TABLE  1 

Procuring  Parameter*  for  Imagery  Example* 


PlfWMtlf 

Halifax  (7m  az)  1-Look 

Trob  RivisrM  (26m  az)  4-Look 

h. 

796.629  km 

796.853  km 

re 

6368.110  km 

6367.623  km 

V«. 

6775.349  ^ 

6774.502  - 
a 

a 

86.629° 

87.860” 

r1 

861.062  km 

866.781  km 

"c 

2s 

1.276  GHZ 

1.276  GHz 

1 

T 

1646.7603  Hz 

1646.7603 

4 

L, 

13.62  km 

4.13  km 

K 

2s 

-0.56276  — 
pare 

-0.56275  jjjgf 

r 

33.9  pmc 

33.9  |t<ac 

*R 

15.7m 

15.7m 

".z 

13£m 

43.9m 

Ar|  (approx.) 

0.26  km 

4.0  km 

V  Vi 

— 

0.516* 

8.  SUMMARY 

A  general  theory  of  two-dlnenalonal  SAR  processing  has  been  presented. 
The  forns  of  the  range  and  aslnuth  radar  signals  returned  fron  a  general 
point  target  have  been  derived,  and  It  has  been  shown  that  these  signals  are 

not  In  general  independent ,  l.e. ,  they  ere  coupled.  An  epproxlaete  closed 
forn  solution  describing  the  inege  produced  by  perforning  e  tvo-dinenslonal 
convolution  of  the  SAR  signal,  returned  iron  s  point  target,  with  the  two- 
dlnenslonal  lapulse  response  of  the  nstchsd  filter  has  also  been  given.  By 
■sans  of  this  solution  it  has  bean  shown  that  a  single  two-dinenslonal 
reference  function  can  be  used  to  produce  high  quality  inages  fron  radar 
slgnala  obtained  over  a  swath  in  range.  However,  these  Inages  oust  be 
geonstrlcaUy  corrected  to  renove  positional  errors  Introduced  by  the 
processing.  Iquatlone  have  been  developed  to  describe  both  the  positional 
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errors  and  the  operations  required  to  remove  them.  Finally,  the  modifications 
required  to  extend  the  theory  to  Include  the  case  of  nulti-look  processing,  or 
noncoherent  averaging  have  been  developed. 
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APPENDIX  A 


i 


i 

! 


Derivation  of  V  „  and  V 
en  ep 

A  brief  derivation  of  (4)  and  (5)  are  given  here.  Refer  to  Figure  1. 
From  rules  of  spherical  trigonometry  for  right  spherical  triangles  one  has 

sin(ir-6 . )cos* 

cos*  ■  — — —  ■  °  .  *  0.  1*  ir/2  rad.  (A.l) 

8  [1  -  sin2Or-0-)sin2*  J*  1 

X  o 

-  0,  6^  “  ir/2  rad. 

and 


sin(ir-0  )  -  sin*,  ./sin*  .  (A.2) 

i  xat  o 


The  velocity  components  are 


V 

en 

and 


V  - 

ep 

where  ue  is  the  planet  rotational  frequency.  After  some  manipulation  and  the 

use  of  sin*  ■  /1-cos2*  ,  (4)  and  (5)  are  obtained, 
s  s 


(«e*e  cos*lat)cos*g  (A. 3) 

(w  r  cos*.  ) sin*  .  (A.4) 

e  e  lat  8 


T 
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APPENDIX  B 


Range  as  a  Function  of  Orbital -Position  In  the  Presence  of  an  Orbital 

Vertical  Velocity  Component 


I 


To  include  a  vertical  velocity  component,  V  in  (17)  h8  is  replaced  by 
h  +  Vvs/V  .  The  second  term  is  the  radial  displacement  from  the  circular 
orbit  undergone  by  the  spacecraft  in  sweeping  through  the  angle  9^.  The  range 
r2  then  becomes 


r0  *  lr  +  Ir  +  h  + 
2  )  e  I  e  s 


(r«  +  h*  +  %*J  -2r«('*  +  h.  +  %°) 


[cos0cos0.  -  sin0cosocos0,  J 

D  D 


(B.l) 


When  the  first  three  terms  of  the  Madaurin  series  are  derived  from  the  equa¬ 
tion  for  r2,  as  in  Section  2,  the  a-coefficients  are  seen  to  be 

(B.2) 


ao‘ri 


sln0  V  V 

a |c  - - £  +  -  tp-  cosOj 

1  rja  tan.  Vfiq  a  Veq 


(B.3) 


2  2r, 


(B.4) 
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APPENDIX  C 


Derivation  of  Equation  Describing  the  Transit  Time 


The  distance  travelled  between  transmission  and 
(see  Figure  11) 


reception  of  a  pulse  is 


(C.l) 


where  rT  and  rR  are  the  ranges  to  the  target  at  the  times  of  pulse  transmis¬ 
sion  and  reception,  respectively.  Then,  following  the  treatment  of  [12],  the 
range  at  the  time  of  transmission  if 


r 


T 


AB 

sin  a 


(C.2) 


where  AB  is  the  distance  between  A  and  B,  and  'a'  is  the  angle  between  B&J. 
and  AS^  as  shown  in  Figure  11.  The  angle  ’b’  between  and  S^A  is  related 

to  ’a'  as  follows 


4 


sin  a  ■  -  cos  b. 


(C.3) 


Therefore,  the  range  at  the  time  of  reception  is  found  from  the 
solution  of 

r»  “  rT  +  I  (r  +r 
R  T  [c  'T  R 

which  can  also  be  written  as 


rR  '  rT  +  [f1  (W]  +  2rT  "I1  <rT+rR)  ^  ‘  <C-» 

If  one  factors  rR-r^  into  (rg~rj) (r  g-Hr^) ,  and  divides  both  sides  by  (rR+r_) , 
the  following  equations  are  obtained  K  T 


>]- 


2rT  (rT+rR)C°8  b 


(C.4) 
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Thus 


V 


+  2-SS. 

C 


At  * 


C 


(C. 


(C. 


Since  (V  /c)  is  always  much  smaller 
the  equat ion  ^ 


than  unity  and  usually  AB  «  r^, 


At  * 


(C 


is  a  very  good  approximation  to  (C.8) 


